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We investigate the electronic band structure and transport properties of periodically alternating 
mono- and bi-layer graphene superlattices (MBLG SLs). In such MBLG SLs, there exists a zero- 
averaged wave vector (zero-fc) gap that is insensitive to the lattice constant. This zero -k gap can be 
controlled by changing both the ratio of the potential widths and the interlayer coupling coefficient 
of the bilayer graphene. We also show that there exist extra Dirac points; the conditions for these 
extra Dirac points are presented analytically. Lastly, we demonstrate that the electronic transport 
properties and the energy gap of the first two bands in MBLG SLs are tunable through adjustment 
of the interlayer coupling and the width ratio of the periodic mono- and bi-layer graphene. 

PACS numbers: 


Since it was first successfully fabricated in experiment 
approximately ten years ago jlj, graphene has become 
an important research topic in condensed matter physics 
and material science. Indeed, because of its unique 
characteristics, graphene has the potential to take the 
place of Si-based semiconductors in future applications 
[ljjlOj. In recent years, many important properties of 
graphene have been explored, such as chiral tunnelling 
[1], the giant carrier mobilities of this material [gj, the 
unusual integer quantum Hall effect I, E3 , and the 
edge-dependant spectra of graphene nanoribbons fTTI . IT^ . 
In particular, for graphene-based superlattices (SLs), 
new Dirac points (DPs) and the zero-averaged wave 
vector (zero-fc) gap have been observed [IHMiTI ] . The 
electronic properties of monolayer graphene (MLG) and 
bilayer graphene (BLG) SLs of various sequences have 
been studied using the transfer matrix method 
In addition, structures with biased potentials 
and heterostructures 24-27| have been investigated. 
However, few studies have discussed the properties of 
periodic mono- and bi-layer graphene (MBLG) SLs, in 
which the MLG is decoupled from the BLG. 



In this paper, we study the electronic band gaps and 
transport properties of MBLG SLs. We demonstrate 
that a zero -k gap and extra DPs exist under certain 
conditions. The zero -k gap opens and closes periodically, 
while the number of extra DPs increases with increased 
lattice constant A values. The relationship between the 
zero -k gap and the interlayer coupling t' B for MBLG 
SLs with periodic sequences is revealed. Furthermore, 
the effects of the average interlayer coupling t', which 
is dependant on t' B and the ratio of the MLG and BLG 
widths (wa and Wb, respectively), on the electronic band 
structures of MBLG SLs are explored. This provides 
a tunable method of controlling electronic conductance 
using MLG and BLG by varying the t' B values (note that, 
for MLG, the interlayer coupling is zero). 




FIG. 1: (Color online) (a) Schematic representation of mono- 
and bi-layer graphene (MBLG) superlattice (SL) periodically 
aligned in x direction. A and B denote monolayer (MLG) and 
bilayer graphene (BLG), respectively. Wa(b) is the width of 
the A (B) region and A = wa + wb is the lattice constant. 
Note that, in the MLG regions, there is no interlayer coupling 
and t'j = 0 in Eq. ©• In the BLG regions, the non-zero f) 
defined in Eq. © occurs, (b) Profiles of periodic potentials 
applied on graphene SL. Va(b) is the square potential applied 
on the A (B) region. 


We begin with the electronic structure of BLG with 
the energy and wave vector close to the K point, such 
that the one- par ticle Hamiltonian for the BLG is given 

by [US,UMll 


/ Vj(x) 7r t'j 0 \ 

7ri Vj(x ) 0 0 

tj 0 Vj( x) 7T t 

V 0 0 7T Vj{x) ) 


(1) 


Here, Vj(x) is the electrostatic potential applied on the 
material, which is constant within each potential barrier 
or well. Further, 7r = — ihvF[-§^ — i-§^\ and 7ri = 
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— Hivf[-§^ + ' L ~§^\ are the momentum operators,where 
vf ~ 10 6 m/s is the Fermi velocity and tj indicates 
the interlayer coupling in the BLG regions. It should 
be noted that tj can be tuned by adjusting the interlayer 
distance [32]. The wave function is expressed by four- 
component pseudospinors $ = (tpi, </> 2 , ^ 3 , Wa) T ■ Note 
that, when tj =0, Eq. © reduces to the MLG case. As 
a result of the translation invariance in the y direction, 
the wave function can be rewritten as <p m = ip m e lkyy , 
m = 1,...,4. By solving this eigenequation, the wave 
functions at any two positions x and x + Ax inside the 
jth potential can be related by a transfer matrix 0 



M 3 = 




and 


M± = 


( cos (qj Ax=p Qj ) 
cos Qj 

. k'j sin(< 7 j Ax) 
^ kj cos Qj 


*f- sin te Aa 0 

cos(qj AxdzQj ) \ 
cos Qj ' 


FIG. 2: (Color online) Electronic band structures for different 

(2) A with wa = wb = (a) 15, (b) 30, and (c) 60 nm, and 
electronic band structures for different potential width ratios 
wa/wb = (d) 1, (e) 3/2, and (f) 2/1 under fixed wb = 15 
nm. The dotted lines denote the centre position of the zero -k 
gap. The other parameters are Va = 100 meV, Vb = 0 meV, 
and t' B = 10 meV. 

( 3 ) 

dispersion at any incident angle follows the relation 


where kj=(E — Vj)/hvF, Ax is the interval of 
any two positions inside the jth potential, q.j = 

sign(kj)^Jkj — ky — tjkj for kj — k 2 —tjkj > 0 , otherwise 

qj = i\J\kj - - tjkj |, tj -)• tj/hv F , kf=k 2 + qj, and 

f lj = arcsin(fc ?/ /fc'). 

If tj = 0, i.e., the BLG is decoupled as MLG, the 
transfer matrix becomes 


M± = 


' cos (qj Atc=p Qj ) 
cos Qj 
^ sinfeAz) 

1 cos Qj 


^ sin(qjAx) 
cos {"2 j 

cos (qj A x±Qj ) 
cos Qj 



( 4 ) 


Compared with the BLG case, for MLG with tj = 0, the 
wave vector kj inside the jth potential reduces to kj . The 
other parameters have the same forms in both cases. 

Using the boundary conditions, the transmission 
coefficient t = t(E, k y ) can be expressed as 


2 cos f 2 o 

(^22e _ * n ° + erne zQ “) — a;i2e i ( ne_n °) — CC 21 ’ 


( 5 ) 


where Xij ( i,j = 1 , 2 ) are elements of the total transfer 
matrix •X’jv==IIjL 1 Afj, with N being the total number 
of potential barriers and wells. The transmissivity is 
T(E,ky) = \t(E,ky)\ 2 . 

Let us consider a periodic sequence of MBLG SLs, 
which is perhaps the simplest sequence (see Fig. [TJ. 
Here, the width and the applied potential of the MLG 
on each unit are wa and Va, respectively, and the 
corresponding BLG parameters on each unit are wb 
and Vb, respectively. According to Bloch’s theorem, 
for an infinite periodic structure (AB) N , the electronic 


cos[/? x A] = - Tr[M A M B \, 


= cos ( q A w A + qsw B ) + sin ( q A w A ) sin ( qBW B ) 

2 k B {q A qB + kj) - k A (k 2 B + k'jj) 


2 ksq A qB 


( 6 ) 


where A = w A +wb ■ Using | cos \/3 x A] | < 1, we can find 
the real solution of j3 x for passing bands. Otherwise, the 
non-existence of real (3 X indicates a band gap. 

Figure[2]demonstrates the properties of electronic band 
gaps under different lattice parameters, indicating that 
there is a zero -k gap in such periodic MBLG SLs. From 
Figs. Ha)-He), it is apparent that the location of the 
zero -k gap is independent of A. In this case, it is 
positioned at approximately 52.63 meV for a variety of 
A. Under the same lattice parameters, the location of 
the zero-fc gap in periodic MBLG SLs is higher and lower 


than those of MLG [l 6 | and BLG [17] SLs, respectively. 
Extra DPs may appear at k y ^ 0 for larger A, as shown in 
Fig. H c ); this will be discussed below. From Figs. H^) - 
Hf), the location of the zero -k gap varies with changes 
in the w a /wb ratio. For example, the zero -k gaps are 
positioned at energy E = 52.63, 62.00, and 68.31 meV 
for w a /wb = 1, 3/2, and 2/1, respectively. 

The location of the zero -k gap is determined by k = 
Y^j=ikj w j/ Y^j=i w j = 3 M- Note that kj should 
be replaced by kj for BLG. When the number of the 
MLG regions is equal to that of the BLG regions, the E 
corresponding to k = 0 can be easily found. Initially, 


E = 


2 (V A w\ - V B w %) - t' B w 2 B 

2K-<4) 


V w a w%(2V a - 2 V B - t ' B ) 2 + t%w 2 B {w B - w A ) 2 


2(w 2 a -w 2 b ) 


(7) 
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Then, for wa = wb, Eq. © reduces to 


E = 


vj-v3- Wb 

2V a - 2V B - t' ' 


( 8 ) 


The locations of the zero -k gaps in Fig. [2] are in 
agreement with Eqs. © and ©. 

We now consider a method to determine the locations 
of the extra DPs, which should obey the equation 


cos \/3 x A] = 1. 


(9) 


Using Eq. ©, when qAWA = — Qb'uJb = m.7r (m is a 
positive integer), Eq. © above is satisfied. For wa = 
wb = w, the zero -k gap will close at the normal incident 
angle {k y = 0) with w satisfying the condition 

mnhv F (2V A -2V B -t' B ) 0 0 , 

v ”" = {v B -v A m-v A + <' B yl m = w - ) (I0) 

At oblique incidences (k y ^ 0) and when qAW a = 
— q B w B = rn7r, the extra DPs appear. They are located 
at 


ky^rn — At 


'(Va-VbYWb-Va + Vb)* 
Vv 2 f ( 2V a ~2 VB-t' B y 


/m 7 T\ 2 
V w ) 


(ii) 

The number of extra DPs can be obtained from Eq. ED 
using the limiting condition: k y < \E/Hvf\. 

Extra DPs exist in the band structures of periodic 
MBLG SLs, which have already been shown in Fig. 
He). In Fig. ©a), the transmission probabilities of 
the centre of the zero -k gap are plotted in order to find 
the extra DPs. In MLG SLs with periodic sequences, 
the charge carriers have perfect transmission at normal 
incidence and a DP is always positioned at the centre 
of the gap 0- However, a different scenario occurs 
in the case of MBLG SLs with periodic sequences. As 
illustrated in Fig. ©a), the zero-fc gap closes with a fixed 
period and extra paired DPs occur as w increases; this 
is characterized by the large transmission probability in 
Fig. ©a). In addition, it is important to notice that 
the DP does not exist at the normal incidence, and that 
it always has large reflection unless the gap is closed. 
From Eq. ED, the gap-closing period is 28.56 nm for 
the parameters given in Fig. ©a). Besides, extra DPs 
do not exist unless A is larger than one period, as can be 
seen in Figs. ©a)-©c). Fig. ©d) shows that the centre 
position of the zero -k gap is independent of A. The effect 
of the BLG’s t' B on the band structure is shown in Fig. 
©e). It is apparent that the width of the zero-fc gap 
increases significantly as t' B increases. The dotted dark 
line denotes the centre position of the zero -k gap and 
satisfies a nonlinear relation [see Eq. (©)]. This differs 
from the linear relation obtained in the case of BLG SLs 
0. 

We next consider the total conductance and Fano 
factor (the ratio between the shot noise power and 



i 


w (nm) 



(d), 

> 

0 

E 


k (nm' 
y v 


851 


80 


S 75 




70" 


20 40 60 

w (nm) 



t' B (meV) 


FIG. 3: (Color online) (a) Transmission probabilities of zero- 
fc gap centre in finite periodic potential structure ( AB ) 25 . 
Electronic band structures for wa = wb = (b) 28.56 and (c) 
57.12 nm. (d) Band-gap structure dependence on w when 
wa = w B = w with k y = 0. (e) Effect of t' B on band-gap 
structure with wa = wb =20 nm. The other parameters are 
Va = 150 meV, Vb = 0 meV, and t' B = 10 meV. 


current) [.‘{.'ll 1351 for different A. In Fig. ©a), the 
conductance of the zero -k gap exhibits different values for 
different A. For fixed A, where wa = = 15 nm, the 

conductance is approximately zero when the zero -k gap 
opens. Specifically, a pair of extra DPs exist at the centre 
position of the zero -k gap for wa = wb =45 nm [see 
Fig. ©a)]. In this case, the conductance of this position 
is slightly higher than zero and the angular-average 
conductance curve forms a linear-like cone around the 
extra DPs. There is no unique DP at the zero -k gap in 
MBLG SLs with periodic sequences, as is the case for 
MLG SL 0,0, 


Therefore, the Fano factor at the 
centre position of the zero -k gap can be removed from 
1/3. 

All the previous results were obtained by focusing 
on the one-mode wave function with longitudinal wave 
vector qj. Experimentally, this scenario can be realized 
using, for example, an adatom to create a wave function 
with a certain longitudinal wave vector. According 
to the result presented in Ref. [36|, the two-mode 


wave function with no mixing can be considered as the 
primary contributor to the conductance or the Fano 
factor [in Ref. [36|, see the values of Tf and TZ and 
compare them with TT and Tf/[. Here, we ignore other 
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FIG. 4: (Color online) (a) Conductance and (b) Fano factor 
vs Fermi energy focusing on one-mode wave function in 
finite periodic sequence (AB) 25 . Averages of (c) conductance 
and (d) Fano factor vs Fermi energy considering two-mode 
wave function in finite periodic sequence (AB) 25 . The other 
parameters are identical to those in Fig. [3] 
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FIG. 5: (Color online) Electronic band structures for 
wa/wb = (a) 2, (b) 1, and (c) 1/2 with fixed wb = 15 
nm. The other parameters are Va = Vb = 0 meV and 
t' B = 20 meV. (d) Energy gap vs t' B . (e) Energy gap vs 
w for wa = wb = w. The double arrows in (a)-(c) denote 
the energy gap widths of the two modes and the dashed lines 
in (d) and (e) denote the relation given in Eq. (1121 . 


complex transport processes that mix different wave 
function modes, and instead consider the contribution 
of the wave function with longitudinal wave vector q' = 

sign (kj)Jkj — ky + t’jkj for — ky + t'jkj > 0, otherwise 

q'j = iyj\k? — k 2 + t'jkj |. The averaged conductance 
and Fano factor produced by both modes of the wave 
functions with qj and q' are plotted in Figs. |U(c) and 
Hd), respectively. In these figures, the conductance 
curves exhibit several dips, whereas the Fano factor 
curves exhibit several small peaks, all of which are 
indicated by red dots. The groups of dips or peaks 
are located at 72.6 and 77.6 meV and are, therefore, 
separated; these separations are primarily caused by 
the zero -k gaps of the wave functions with q' and qj. 
For the case of wa = wb = 15 nm, the average of 
the conductance at approximately 77.6 meV is increased 
compared with the zero conductance in Fig. [2a), and the 
position of the conductance minimum is changed from 
77.6 to 75.0 meV. Based on the above analysis, it may 
be concluded that the result obtained by focusing on the 
one-mode wave function is meaningful as a reflection of 
the transport properties of normal scenarios involving the 
two-mode wave function. 

Finally, we consider another interesting effect. In the 


absence of an applied potential on the periodic MBLG 
SLs, the energy gap of the first two bands can be regularly 
tuned by controlling the wa/wb ratio. The yellow 
regions in Figs. Eta)-Etc) are occupied by electrons of the 
propagating mode, whereas the dark regions are occupied 
by electrons of both the propagating and evanescent 
modes. As illustrated in these figures, the gap width 
of the first two bands changes from tfl/2, to 2t' B /2> 

for wa/wb = 2, 1, and 1/2, respectively, under the given 
parameters. Further, the gap width is proportional to 
the wb/{wa + wb) ratio and roughly obeys 


t' B WB 

wa + w b ' 


( 12 ) 


Figs. Md) and [5[e) show the correct range of Eq. 
da. From Fig. [5f d), it is apparent that Eq. (fT2l) 
accurately describes the regularity when t' B is smaller 
than approximately 50 meV in the case of w B = 15 
nm, and the deviation is larger for larger t' B . A also 
affects the correctness of the above relation, as shown in 
Fig. EJe). It is concluded that the actual gap width is 
smaller than the approximated result given by Eq. m, 
which indicates that, as t' B or w increases, the effect of 
the MLG becomes stronger than that of the BLG. Based 
on the similarity between the band structures in Fig. [5] 
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and those of the BLG, we may define E g as t'. As t' 
increases, the band structure becomes more similar to the 
parabolic band structure of BLG. For a periodic structure 
with alternating interlayer couplings in total two-layer 
graphene, Eq. IT2l) (which concerns the gap width) can 
be modified to 


t' A w A + t' B w B 

^9 — i 

WA + W B 


(13) 


where t' A is the interlayer coupling in the A regions, 
which are also BLG. The t' in MBLG SLs corresponds 
to the border between the propagating and evanescent 
wave modes, and this BLG border has been shown in 
Duppen and Peeters’ interesting work (see Ref. i36). 
Furthermore, this border is tunable through adjustment 
of t' A and t' B and wa/wb- Therefore, this may constitute 
a new method of controlling the electronic conductance. 

In summary, we have studied the electronic band 
structures and transport properties of mono- and bi¬ 
layer graphene (MBLG) superlattices (SLs) with periodic 
sequences. The zero -k gap is robust against the lattice 
constant A = wa + w B , but sensitive to both the ratio of 
the potential widths wa/w B and the interlayer coupling 
t' B . The analytical condition determining the locations 
of extra Dirac points (DPs) ( k y ^ 0) was presented 
and it was shown that the number of DPs increases as 
A increases. Moreover, the conductance of the charge 
carriers in periodic MBLG SLs was investigated, and it 
was shown that the angular-average conductance curve 
forms a linear-like cone around the extra DPs, similar to 
that observed in the case of MLG SLs M- If no potential 
barriers are applied, the energy gap width of the first 
two bands or the average interlayer coupling t' is tunable 
through adjustment of t' B and wa/wb- This finding is 
important as regards research on interlayer coupling in 
two-layer graphene, and may constitute a new method of 
controlling the electronic conductance. Further related 
results will be reported in the near future. The findings 
of this and future studies will be useful for future research 
and applications concerning MBLG SLs. 
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